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ABSTRACT 

The star formation rate (SFR) and black hole accretion rate (BHAR) functions are measured to be propor- 
tional to each other at z < 3. This close correspondence between SF and BHA would naturally yield a BH 
mass-galaxy mass correlation, whereas a BH mass-bulge mass correlation is observed. To explore this ap- 
parent contradiction we study the SF in spheroid-dominated galaxies between z = 1 and the present day. We 
use 903 galaxies from the COMBO-17 survey with M, > 2x 1O 1O M , ultraviolet and infrared-derived SFRs 
from Spitzer and GALEX, and morphologies from GEMS HST/ACS imaging. Using stacking techniques, 
we find that <25% of all SF occurs in spheroid-dominated galaxies (Sersic index n >2.5), while the BHAR 
that we would expect if the global scalings held is three times higher. This rules out the simplest picture of 
co-evolution, in which SF and BHA trace each other at all times. These results could be explained if SF and 
BHA occur in the same events, but offset in time, for example at different stages of a merger event. However, 
one would then expect to see the corresponding star formation activity in early-stage mergers, in conflict with 
observations. We conclude that the major episodes of SF and BHA occur in different events, with the bulk of 
SF happening in isolated disks and most BHA occurring in major mergers. The apparent global co-evolution 
results from the regulation of the BH growth by the potential well of the galactic spheroid, which includes a 
major contribution from disrupted disk stars. 

Subject headings: galaxies: evolution — galaxies: active - quasars: general 



1. INTRODUCTION 

The last decade has seen the discovery and character- 
ization of an unexpectedly tight correlation between the 
mass of supermassive black holes (SMBH; Mbh) and the 
mass ( MbnWe) or velocity disp er sion of their host galaxy's 
bulge (Magorr ian"etal.l 119981: iFerrarese & Merrittl 120001: 
Gebh ardt et all l2000t iMarconi & Hunt! 120031: lHaring & Rixl 
2004). The black hole mass appears to be mo st strongly corre- 



lated with the bulge mass, not the total mass jKormendvl200U 
lHoll2007l) : the scatter in this relation between bulge mass and 
black hole mass is e stimated to be less than a factor of two 
(lHaring & Rixll2004h . 

This relationship indicates that galaxy and black hole for- 
mation and evolution are interconnected. 8 In its weak form, 
such an interconnection could result if the black hole growth 
is limited by the wider galaxy environment ('dog wagging 
the tail'). Stronger forms of interconnection are also possi- 
ble. The energy released by a growing SMBH is sufficient, 
if it couples effectively with its surroundings, to have dra- 
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M bulge and M B h- 



matic consequences on the galaxy ('tail wagging the dog'). 
For example, the injection of energy into hot halo gas in 
galaxy clusters is indicated by the inflation of bubbles a nd the 
propagation of sound waves dMcNamara et alj l2000l 12001 
iFabianet all 120031: iRaffertv et alj 120061) . Furthermore, it is 
possible that SMBHs with high accretion rates drive power- 
ful outflows that remove (at least some) cold gas from galax- 
ies (IChartas et alJl2003t ICrenshaw et al.ll2003t iPounds et alJ 
120031: iTremonti et alJl2007fT In fact, feedback from an accret- 
ing SMBH/AGN may help address many of the most diffi- 
cult issues affecting models of galaxy evolution in a cosmo- 
logical contest: e.g., overcooling in massive halos and late 
star formation in elliptical galax ies (Binnev & Tabor! [1995 : 
I Croton et all 12006b Bower et all 120061: iMonaco etafl 120071: 
ISomerville et alJl2008h ~ 

A diversity of models has been constructed to explore the 
interrelationship between bulge mass and black hole mass, 
featuring bot h weak and strong aspe cts of the possible inter- 
connections (ISomerville et alJfeOOSl, and references therein). 
However, the details of the physical processes governing BH 
formation and growth, and the physical mechanisms whereby 
SMBHs couple with their host galaxies and surroundings re- 
main poorly understood. It will be impossible to model the 
full range of interconnected processes from theoretical first 
principles for some time, as the physics of the formation of 
SMBHs and galaxy evolution spans at least a 7 orders of mag- 
nitude in scale. Empirical and phenomenological constraints 
are therefore extremely important. 

A common feature of most "unified" models of SMBH and 
galaxy formation is the concept of co-evolution, whereby the 
SMBHs and their galaxies evolve under each other's influ- 
ence. The term co-evolution means many things to many 
people (e.gJSilk & Reesll998tlKa"u ffmann & Haehneltl200oi: 
I Wvilhe & Loebl I2003t iGranatoetalJ l2004t iMerloni et alj 
20041: lHaiman Ciotti & Ostrikeril2004t Icattaneo et all 120051: 
Fontanot et al. 2006). For the purpose of having strawmen to 
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take aim at, we outline three possible scenarios: 

• Strong Co-evolution: In its strongest sense, co- 
evolution implies that bulges and black holes grow to- 
gether, i.e., in the same objects and at the same time. 
In this case, BH growth could be thought of as a sort 
of a 'tax': whenever star formation occurs, the SMBH 
accretes some fraction of the mass involved. This sce- 
nario in its simplest form can already be ruled out — 
there are many examples of star-forming bulges with- 
out AGN activity and bulges with AGN activity and no 
SF (e.g., M87). 

• Time Offset: One could imagine a scenario in which 
bulge stars and SMBHs are formed in the same event, 
but on different timescales and/or stages of the event. 
For example, star formation might occur predominantly 
in the early stages of a merger, while BH growth might 
occur in the later stages (e.g., iDi Matteo et alj 120051: 
Hopk ins et al.l l2005, 2006). Such a picture would pre- 
dict little evolution in the BH-bulge mass correlation 
for an unbiased ensemble, but galaxies in the throes of 
their bulge building events (or black holes in their most 
rapid growth phase) could exhibit significant deviations 
from the scaling relations; this hypothesis is consider- 
ably more challenging to rule out. 

• Regulated Growth: A third possibility is that the 
major episodes of SF and BH accretion occur in dif- 
ferent events, but the growth of one component reg- 
ulates that of the other. For example, if black hole 
growth is regulated by feedback fr om the AGN itself , 
as suggested by many workers (e.g..[Silk & Reeslll998t 
Murr ay. Ouataert & Thompson] 120051: [Di Matte o et alJ 
2005), the potential well of the pre-existing bulge de- 
termines how much energy is needed to stop further ac- 
cretion and halt the growth of the black hole, and hence 
the final black hole mass. It is possible, even likely, 
that the bulge stars that dominate this potential well are 
formed mostly in previous star formation "events", i.e. 
in progenitor disks which have subsequently merged. 
In this picture, one would expect a larger possible dis- 
connect between observable episodes of SF vs. BH ac- 
cretion activity, even while the endpoints of this activity 
(bulges and black holes) are expected to end up being 
tightly related. 

There are a number of ways to attack the problem of co- 
evolution observationally. Obviously, measuring the BH- 
bulge mass correlation at the present day at extremely 
low or high mass will provide powerful insig hts into the 
proce s ses linking galaxy and SMBH evolution dBarth et ail 
12005b lLauer et al.M2007al) . The study of the BH-bulge 
mass correlation during BH or bulge-buildin g events is 



times than those seen today (|Treu et alj l2004t IWoo et all 
20061 iPeng et al] l2006t iShields et all [2006. see also, e.g.. 



Alexander et al. l2005f for the smaller ratios of black hole to 



also an important avenue of approach (e.g., Barth et al 



20051 : iGreene & Hdl2006t iTreu et alj|2004 12007b IWoo et al 
20061: IPeng et al.l l2006t iBorvs et all l2005l)~ The most 



commonly-used technique in the latter approach is study- 
ing the masses of actively-accreting Type I AGN for 
which BH masses can be estimated using an empirically- 
calibrated combination of luminosity and broad line region 
line width dWandel. Peterson & Malkanll999tlBaskin & Laoil 
l2005t iKaspi et all 12005b IVestergaard & Petersonl 120061 and 
references therein); current results are controversial but ten- 
tatively support larger BH to bulge mass ratios at early 



bulge mass in rapidly-star-forming galaxies detected in the 
submm). 

With current generations of wide and deep cosmolog- 
ical surveys, another means has opened up: compari- 
son of evolution of the BH accretion rate (BHAR) and 
that of the SFR of ga laxies ([Si lverman et al. 200a see 
iBoyle & Terleviclil [1998 and Franceschini et alJ 119991 for 
early efforts). The history of BH accretion has been ex- 
plored using primarily optical and X-ray techniques, with im- 
portant insights being gained from observations i n the radio 
and i n frare d (Mivaii, Hasinger & Schmidt 2000; [Ueda et al 
lHas inger. Mivaii & Schmid t] 12005b lLa Franca et al 
Bareer et alJ 12005b iDonlev et all 12005b iBrown et al 
iRichards et al.l l2006ah . Furthermore, a decade 
intensive multi-wavelength and redshift surveying has 
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placed powerful and interesting constra ints on the volume 
avera g ed SF history of the Unive r se dHopkins & Beaco~ 
2006: iPerez-Go nzalez et~aT] 120051: iLe Floc'h et al] 120051 
Schim inovich et all l2005). In both cases, the integrated BH 
accretion and SF histories give values that match to within 
a factor of 2 or 3 of the present-day BH mass and stellar 
mass 9 densities dYu & Tremamd 120021 iMarconi et al] 12004b 
iBorch et al.ll2006l: iFardal et al.ll2007blWilkins et al.ll2008l) : we 
will discuss in detail uncertainties and possible inconsisten- 
cies between the integral of the SFR and the present-day stel- 
lar mass in $2] 

The observational maturity of current datasets at z < 1 al- 
lows the use of a novel angle of attack on the problem of co- 
evolution. Recent multi-wavelength/HST surveys have cov- 
ered enough volume that SFR density and BHAR density can 
be estimated in well-defined subsamples of galaxies for the 
first time, offering potentially decisive insight into not only 
the statistical evolution of these quantities, but into whether 
black holes and bulges are actually growing in the same ob- 
jects and at the same time. This is the approach we will adopt 
in this paper. To do this, we first examine the global SFR 
and BHAR (©, then the multiplicity functions 10 of SFR and 
BHAR as a function of redshift (@. In $H we present new 
measurements of SF activity in mass and morphology-limited 
samples of galaxies, and examine accretion activity driven by 
massive BHs in massive spheroids under assumptions widely 
used in literature. In $5] we discuss the results in the context 
of co-evolution and summarize our conclusions. Throughout 
the paper we assume a cosmology with Ho = 70 km s _1 Mpc~' , 
fi M = 0.3 and Q A = 0.7. 

2. THE COSMIC STAR FORMATION AND BH ACCRETION 
HISTORIES 

We first attempt to give some context to the analysis that 
will follow by exploring the cosmic star formation and black 
hole accretion histories. 

2.1. The cosmic star formation history 

During the last decade, a number of studies have con- 
tributed to determining the average SFR per co-moving vol- 

9 In the case of stellar mass density, recycling of nearly 1/2 of the initial 
stellar mass back into interstellar gas must be accounted for to avoid dramat- 
ically overproducing the present-day stellar mass density. 

10 The comoving number density of sources with a given SFR or BHAR 
— e.g., a luminosity function is a multiplicity function 
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Fig. 1 . — Top: Comparison between the volume-averaged SF history and the volume-averaged BH accretion his tory; Bottom: The same comparison but only 
in the redshift range < z < 1. The gray da t a poin t s come from the com pilati on of availab l e mea surements by Hopkins & Beacom 12006). The black squares 
are radio m easurements from Seymour et al. 12008), Smolcic et al. (2009) and Dunne et al. (2008). The black diamonds represent IR+UV measurements from 
IBell et at] 120071) . All measurements are converted to a Kroupa IMF and the same cosmology adopted here. The BH accretion history is derived from AGN 
bolometric luminosity functions given in HRH07 by assuming a radiation efficiency e = 0. 1 . The BH accretion history is shifted upwards by 3.3 dex (a factor of 
2000). 
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ume at different cosmic epochs with various SFR estima- 
tors (e.g., UV, Ha, IR and radio). Th e available measure- 
ments from the literature compiled by Hopkin s~& Beacornl 
(2006), and co mplem ented with recent measurements from 
Seym our et all (12008b . iSmolcic et all (120091) and Dunne et alj 
(2008) are shown in Figure Q] We include also recent esti- 
mates from the use of UV and 24pm data to acco unt for both 
unobscured and obscu red star formation at z < 1 (Zheng et al. 
120071; iBell et al.ll2007h . Here all measurements are collected 
to our adopted cosmology and a Kroupa (2001) initial mass 
function (IMF). It can be seen from FigureQ]that the cosmic 
SFR density increases rapidly with redshift, peaks at z ~ 1.5, 
then becomes flat or declines at z > 1.5, although the un- 
certainties are large at these early epochs. The cosmic SFR 
density decreases by an order of magnitude from the z = 1 to 
the present day, following psfr ~ (1 + z) 2 - 9±02 . 

A complementary measurement to the cosmic SFR density 
is the cosmic stellar mass density. The assembly of stel- 
lar mass should be consistent with the integral of the cos- 
mic SFR density. Current measurements show a reasonably 
good agree ment between the two, at least at redshifts less than 
about unity dBorch et al.H2006t IBell et al.ll2007t iWffldns et all 
120081): about 40-50% of local stars were formed since z - 1 
dDickinson et al.ll2003l:lFontana et al.l2003l:lDrorv et al.H2005t 
iBorch et al.l 120061 ; iRudnick et al. 2006). As an aside, it is 
worth noting that there is a possible mismatch between the 
integral of the cosmic SFR density at z > 1 with the stel- 
lar mass observed to be in place at z = 1: it appears that 
the integral is a factor of ~ 3 in excess of the stellar mass 
formed by z = 1 . The origin of this mismatch is currently not 
well-understood. It is possible that this discrepancy signals 
a break-down in the utility of SFR indicators or a change or 
break-down of a universally-applicable stellar IMF. It should 
also be noted that at z > 1 the estimates of SFR (and stellar 
mass) are highly uncertain; in particular, the relatively sensi- 
tive UV measurements need to be corrected for substantial 
(factors of a f ew) dust extinction using defensible but un- 
certain recipes (Reddv et al. 2008), and measurements of the 
rest-frame IR a nd radio probe on ly the most luminous sys- 
tems (see, e.g., Dun ne et al.ll200 8 for substantial progress to- 
wards this goal using stacking of radio data; it is interesting 
that their resulting SFR densities are much lower than those 
reported previously at z > 1.5). Yet, for our purposes, the pos- 
sibility of a mismatch between SFR and stellar mass at z > 1 
is of relatively little importance; the focus of our work is at 
the better-constrained z < 1 redshift range. 

2.2. The cosmic BH accretion history 

Observed (optically "bright") BH accretion apparently can 
account for nearly all of the BH mass seen in remnants 
today, indicating that the bulk of BH mass growth occurs 
during a lumin ous AGN phase with a radiation efficiency 
e ~ 0.1 (e.g lYu & Tremaind 120021; IMarconi et all l2004t 
iShankar et al j 120041) ! The observed AGN luminosity func- 
tion is therefore a reasonable probe of the cosmic BH ac- 
cretion history. ^From deep cosmological surveys performed 
with modern observational facilities, AGN luminosity func- 
tions have been determined out to z ~ 4 i n hard X-ray 
(e.g. lUeda et al.l 120031; lLa Franca etail 120051: iBarger et al.l 
20051). soft X-ray (e.g.. iMiyaii. Hasinger & Schmidt I I2000T 
Hasinger, Mivaii & Schmidt 2005), optical (e.g.. I Croom et alj 
20041: iRichards et al.ll2006bt[Fontanot et al.120071) an d mid-IR 
bands (e.g jBrown et al.l2 006; Richards et al. 2006a). The in- 
trinsic spectral energy distribution (SED) of AGN is thought 



to be only a function of luminosity (e.g., Mar coni et alj 2004); 
and accounting for obscuration, the AGN luminosity func- 
tions of diffe r ent b ands are correlated through the SED. 
iHopkins et al.l (l2007bl hereafter HRH07) presented bolomet- 
ric AGN LFs derived from the combination of AGN LFs mea- 
sured in multi-wavelength bands, from mid-IR through hard 
X-ray (see their paper and references therein for details of 
the AGN LFs and related uncertainties). We caution that 
Compton-thick AGN are not counted in the HRH07 bolomet- 
ric LFs. Correction for the missing accretion is suggested to 
be < 0.15 dex (a factor of 1.4; see HRH07 for more discus- 
sion). HRH07 presented several evolution models which fit 
the AGN bolometric luminosity function at different cosmic 
epochs. We adopted the best-fit pure luminosity evolution 
(PLE) model and the luminosity-dependent density evolution 
(LDDE) model of the bolometric LF <£(L bo i,z) from HRH07 
to calculate the BHAR per co-moving volume as follows: 

Pbh(z)= f (1 " c) , Lfaol a>(L bo i,z)dL fao i, (1) 
Jo « 2 

where Lboi is the bolometric luminosity and we have adopted 
e = 0.1. By default Lboi is given in units of ergs -1 . We inte- 
grate the bolometric luminosity over the range 43 < log Lboi < 
49. The cosmic BHAR density peaks at z ~ 1.9, and dra- 
matically decreases to the present day. At z > 1 -9, the cos- 
mic BHAR density is likely to decline or stay flat to earlier 
cosmic times altho ugh the uncertainties are large (see also 
Barger et af] 120051) . Roughly speaking, the cosmic BHAR 
density increases as p B H ~ (1 + zf and a ~ 3.4 up to z = 1. 
Considering the uncertainties of different band LFs and un- 
certainties in bolometric correction, a typical error of ~0.6 is 
adopted for the power index a. 

2.3. Comparison between the cosmic star formation and BH 
accretion histories 

Figure[T]shows the cosmic SFR density to be in good agree- 
ment with the cosmic BHAR density scaled by a factor of 
2000 over a wide redshift range from z ~ to z ~ 5. The 
factor 2000 is the ratio of the integral of the cosmic SFR den- 
sity to the integral of the cosmic BHAR density within the 
redshift range from z — to z — 1, where both quantities are 
well determined. In the redshift range < z < 1 a quantitative 
comparison between the cosmic SFR density and BHAR den- 
sity can be made: they increase with redshift with almost the 
same slope of p ~ (1+z) 3 within the errors. Figure[T]suggests 
that the cosmic BH accretion history is parallel to the cosmic 
SF history. 

If the strong co-evolution hypothesis is correct, the shift be- 
tween the volume-averaged SFR and volume-averaged BHAR 
(a factor of 2000; 3.3 dex) should be consistent with the 
difference between the local stellar mass density and BH 
mass density. The integral of the SFR needs to be cor- 
rected for the recycling of mass back into gas during the pro- 
cess of stellar evolution; with our adopted IMF, almost 50% 
of the initially-formed stellar mass is returned to the ISM 
within 5-7 Gyr. Thus, one expects a ratio of local stellar 
mass density to BH mass density of ~ 1000. This is con- 
sistent with the measurements of l ocal stellar mass density 
p.,o = 3.0^* x 10 8 M o Mpc- 3 (e.g.. lBelletani2003l) and lo- 
cal BH mass densit y p BH o = 4.6+J^ x 10 5 M Q Mpc- 3 (e.g., 
IMarconi et al1l2004h . 

Thus, somewhat remarkably, it appears that the histories of 
SMBH accretion and star formation are very similar in shape, 
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TABLE 1 
Best-fit parameters for the SFR 

FUNCTIONS." 



-'■start 




log SFR* 
(Moyr 1 ) 


log<D* 
(MpcT 3 dex" 1 ) 


0.2 


0.4 


1.09 


-2.80 


0.4 


0.6 


1.11 


-2.77 


0.6 


0.8 


1.12 


-2.63 


0.8 


1.0 


1.15 


-2.90 



a The SFR functions are given in the form 
of "double power-law", logd> = log®* + 
a \og S FR, with a = -0.6 and a = -2.2 for 
S FR below and above S FR* , respectively. 



offset by a ratio consistent with the ratio between BH mass 
and stellar mass density at the present day. Put differently, 
it appears that the ratio of stellar mass to SMBH mass den- 
sity is independent of redshift, in seeming accord with the 
co-evolution picture. 

3. STATISTICAL LINKS BETWEEN SF AND BH ACCRETION 
EVENTS 

The overall BH accretion/SF mass ratio results from the 
sum of individual BH accretion/SF events in the same cos- 
mic epoch. Here, a BH accretion "event" means a SMBH in 
the active phase (i.e., AGN) and a SF event refers to a galaxy 
of a given SFR. In this section we address whether the inten- 
sity of star formation events is statistically correlated with the 
intensity of BH accretion events. 

Following the description in 92.21 we convert the AGN 
bolometric luminosity functions of HRH07 into BHAR func- 
tions by assuming e = 0.1. We take the BHAR functions 
described by a double po wer-law as des cribing the statistics 
of BH accretion events. IB ell et all (120071) estimated the SFR 
using bolometric (UV+IR) luminosity for a sample of 7506 
galaxies and derived SFR functions in four redshift slices be- 
tween z - 0.2 and z = 1- Motivated by local IR lumi- 
nosity functions that are w ell-fit with a "double power-law" 
shape (Sa nders et al.ll2.003l having a form of O(L) oc L a with 
a = -0.6(+0.1) and a = -2.2(+0.1) forL < L* and L > L* 
respectively), we re- fit the SFR functions of lBell et alJ (120071) 
with a double power-law. The best-fit parameters for these 
SFR functions are given in Tablefl] We adopt these SFR mul- 
tiplicity functions to describe the statistics of SF events. 

Figure|2]shows the adopt ed B HAR functions and SFR func- 
tions. We have shown in 92.31 that the cosmic BH accretion 
history tracks the cosmic SF history through a universal scal- 
ing factor of 2000. Comparing the BHAR function with the 
SFR function, we find that the former always tracks the latter 
in all four redshift slices from z = 0.2 to z = 1 after being 
re-scaled by a factor of 20 in BHAR and by 100 in num- 
ber density, as shown in Figure [2] The two scaling factors 
20 and 100 are empirically determined to make the agree- 
ment between the two functions as good as possible, partic- 
ularly in the wide regime around the SFR function "knee" 
(i.e., 0< log{SFR/ M yr~') <2), where the vast majority 
of the total SF occurs. We split SF events by their inten- 
sity into three classes: high-intensity (S FR > 100 M yr~'; 
i.e., ultraluminous IR galaxies), medium-intensity (10 < 
SFR < 100 M yr _1 ; i.e., luminous IR galaxies) and low- 
intensity (SFR < 10 M yr _1 ; i.e., "normal" galaxies). Sim- 
ilarly, we split BH accretion events into the same three 
classes by replacing SFR with 2QxBHAR. The three classes 



of AGN roughly correspond to luminous quasars, quasars 
and Seyferts/low-luminosity AGN, respectively. From the 
adopted SFR/BHAR functions, we calculate the volume- 
averaged SFRs/BHARs contributed by the three classes, re- 
spectively. The results are shown in Figure [3] 

As can be seen from Figure [3] the volume-averaged 
SFR/BHAR contained in high-intensity star formation or BH 
accretion events decreases dramatically from z = 1 to z = 0.2, 
whereas the decrease of the volume-averaged SFR/BHAR 
with decreasing redshift is gradually slower for the medium- 
intensity and low-intensity SF and BH accretion events. The 
low-intensity SF and BH accretion events dominate the cos- 
mic SFR/BHAR at low-z and the high-intensity star for- 
mation/BH accretion events start to dominate the cosmic 
SFR/BHAR at z > 0.8. It is clear that for either the medium- 
intensity or the low-intensity class, the volume-averaged SFR 
matches the volume-averaged BHAR remarkably well over 
the redshift range 0.2 < z < 1. For the high-intensity 
class, the agreement is poor, possibly because of the large 
uncertainties in the SFR and BHAR functions in this range 
(SFR > 100 Moyr 1 or 20 x BHAR > 100 Meyr 1 ). We 
conclude that SF events statistically track BH accretion events 
over 0.2 < z < 1 when split by intensity. 



4. CHARACTERIZING THE HOST GALAXIES OF STAR 
FORMATION AND BH ACCRETION 

In §2 and §3, we presented evidence that both the distribu- 
tion (multiplicity functions) and integrals of the SF rate and 
supermassive black hole accretion rate evolve similarly at all 
redshifts, with an offset in integrated rate of a factor of ~ 2000 
at all z < 1. Taken at face value, such a coincidence may 
lead one to predict a correlation between galaxy mass and 
supermassive black hole mass, offset by a factor of ~ 1000 
in zero point (where the expected offset is a factor of 1000 
rather than 2000 because of the recycling of stellar mass dur- 
ing the course of stellar evolution). Such a correlation is not 
observed; instead, it appears t hat bulge mass an d supermas- 
sive black hole mass correlate (Kormendv 200ll). 

In this section, we present new measurements of the re- 
lationship between SFR and BHAR for spheroid-dominated 
galaxies, in order to investigate the seeming disconnect be- 
tween a BH-bulge mass relation on one hand, and the close 
correspondence between SFR and BHAR statistics on the 
other. 

4.1. Links between SFR, stellar mass and morphology 
4.1.1. The data 

We address this issue using a deep, wide-area multi- 
wavelength dataset from the extended Chandra Deep Field 
South (ECDFS). We use optical photometry, photo metric red- 
shift catalogs ( fe/(l +z) ~ 0.02;lWolf et al.l l2004) and stellar 
mass estimates (Bor ch et al.ll2006l) from the COMBO- 17 sur- 
vey for ~ 9000 galaxies with aperture magnitudes wr < 23.5 
mag and z < 1.1 in the 30'.5 x 30' ECDFS. High-resolution 
((X'07) HST imaging from the Galax y Evolution fro m Mor- 
phology and SEDs (GEMS) survey dRix et al.ll200H) covers 
around 800 square arcminutes of the ECDFS in the F606W 
and F850LP passbands, providing optical morphologies for 
~8000 galaxies. A two-dimensional light distribution analysis 
is performed on the F850LP imaging data using the software 
tool GALFIT, providing Sersic index n as a broad measure of 
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the concentration of a galaxy (HauBler et al. 2007|). 11 

Deep far-ultraviolet (FUV; 1350-1750A) and near- 
ultraviolet (NUV; 1750-2800 A) images centered on the 
ECDFS were obtained by the Galaxy Evolution Explorer 
(GALEX: Martin et al. 2005), with a field of view of a square 
degree. The FUV and NUV images have a typical resolution 
~ 5", and a depth of 3.63 /iJy at the 5 <x dete ction level. Data 
reduct ion and source detection is described in Mo rrissev et al.l 
(2005). Deep 24yum data of th e ECDFS were ta ken as part 
of the MIPS GTO observations dRieke et al. 2004). A mosaic 
image was produced with a rectangular field of 90' x 30'. 
The 24 fim image has a point spread function (PSF) with full 
width at half maximum (FWHM) =* 6". Sources are detec ted 
down to 83 /iJy at the 5 cr level (see lPapovich et af]|2004 for 
details of data reduction, source detection and photometry). 

4.1.2. Sample selection 

The goal of §4 is to study the SFR in massive galax- 
ies split by morphology (or, more precisely, structure, as 
we use Sersic index n to differentiate between spheroid- 
and disk-dominated galaxies). Accordingly, in this section, 
we define the mass-limited, SFR-limited and structurally 
(morphologically)-selected galaxy samples that we use for 
further study. 

A sample of 903 massive galaxies (M, > 2 x 10 10 M G ) in 
the redshift range of 0.2< z <1 is selected from the COMBO- 
17/GEMS survey in the ECDFS. These objects are located in 
the overlap area (~600 square arcminutes) of the HST, GALEX 
and MIPS observations that allow for SFR measurements. We 
have removed X-ra y point sources dete cted in the Chandra 
250 ks observation (Leh mer et al.l 12005) in order to exclude 
AGN-heated dust and AGN UV emission (we describe later 
how this selection criterion is of only modest importance; less 
than 15% of star formation is in such systems). For the present 
comparison, we identify spheroid-dominated galaxies with 
galaxies of Sersic index n > 2.5. The criterion n - 2.5 is a 
reasonable and reproducible separation betw een concentrated 
early-type galaxies a nd late-ty pe galaxies dBell et aljf 2004b; 
iPannellaet all 120061 see, e.g.. Ivan der Wei et al.l 120081 for a 
discussion of the relationship between morphology and Sersic 
index in the SDSS). In what follows, we use this 903-galaxy 
sample to define three galaxy subsamples: 1) mass limited 
(M» > 2 x 10 10 M o ; 903 galaxies); 2) mass and SFR limited 
(SFR > 3M yr _1 ; 308 galaxies, SFRs are defined in §4.1.3); 
and 3) mass and morphology-limited (Sersic index n > 2.5; 
493 galaxies). 

4.1.3. Estimating SFR from the UV and IR 

We divide the galaxies of each subsample into four even 
(Az = 0.2) redshift slices. Average UV and IR luminosities 
are estimated for galaxies in each redshift bin and each sam- 
ple. We use stacking techniques to estimate the average fluxes 
of individually un detected objects in a subset of galaxies (see 
Zheng et al. 2006 for details about the 24 /mi stacking). This 
is important for obtaining a complete SFR estimate for mas- 
sive spheroid-dominated galaxies because most such galax- 
ies are intrinsically faint in the UV and mid-IR. The stacked 

1 1 Surface brightness dimming is likely to introduce a redshift-dependent 
systematic effect. Disk components (typically of relatively lower-surface 
brightness) will probably be less frequently detected at higher redshifts (see 
IShi et alj|2009t for a related discussion). A tendency towards a higher frac- 
tion of disk contamination in higher-redshift spheroids could result, affecting 
the inferred redshift evolution of the SFR density in n > 2.5 systems. 



TABLE 2 

Average SFR and Stellar mass for three galaxy sub-populations 
over < z < 1 ■ 



z N obj < SFR > <M,> psFR 

(M yr-') (10 l0 M o ) (1(T 2 M Q yr" 1 Mpc~ 3 ) 



massive galaxies 



0.04 








0.4±0.2 


0.3 


64 


3.1±0.2 


5.0±0.5 


1 . 1 +0.5 


0.5 


230 


4.3±0.2 


4.5±0.6 


1.2±0.4 


0.7 


421 


8.3±0.3 


5.7±0.5 


2.5±0.7 


0.9 


188 


15.1+0.8 


5.5±0.7 


3.5+1.9 






massive spheroid-dominated 


galaxies 


0.04 










0.3 


35 


0.9±0.1 


4.5±0.8 


0.22±0.12 


0.5 


122 


1.3+0.1 


6.4±0.8 


0.22±0.07 


0.7 


249 


4.2±0.2 


6.7±0.7 


0.65±0.19 


0.9 


87 


9.1±0.8 


6.8+1.2 


1.01±0.59 






massive S FR > 3 M yr 1 ; 


»alaxies 


0.04 








0.2±0.1 


0.3 


19 


7.1 


4.1 


0.7±0.4 


0.5 


89 


10.2 


4.9 


1.1 ±0.6 


0.7 a 


176 


15.2 


4.5 


>1.9±0.8 


0.9 a 


81 


23.0 


4.6 


>2.3±1.4 






massive SFR> 10M o yr' 


galaxies 


0.04 








0.04±0.03 


0.3 


3 


14.3 


5.9 


0.2±0.1 


0.5 


29 


20.2 


4.9 


0.7±0.2 


0.7 


99 


22.0 


4.8 


1.6±0.5 


0.9 


63 


27.4 


4.0 


2.1 + 1.2 



11 Sample is incomplete in this redshift bin. The average SFR is 
overestimated and the volume-averaged SFR is underestimated. 

fluxes and individually detected fluxes are combined to obtain 
the mean flux for the given subset of galaxies. This procedure 
was applied to MIPS 24 /jm and GALEX FUV and NUV. Aver- 
age FUV, NUV and 24 fim luminosities are obtained for each 
subset of galaxies. Optical photometry in the U, B, V,R and 
/ bands from the COMBO- 17 survey is available for all in- 
dividual sample galaxies, allowing calculation of the average 
luminosities in these bands. 

The total UV luminosity is estimated by integrating the 
spectral energy distributions over the wavelength range of 
rest-frame 1500 - 2800 A from linear interpolation of the 
FUV, NUV and optical U,B, V,R and / bands. We estimate 
the total IR luminosity from the observed 24 yizm luminosities 
using three set s of luminosity-dependent IR SEP templates 
(Laga che et all l2004t iDale & Heloul l200l iCharv & Elbazl 
2001). 12 Although these templates are derived from lo- 
cal star-forming galaxies, they can be used to represent IR 
SEDs for distant s t ar-forming galaxies (|Marcillac et alll2006t 
Zhen g et alj [20071: iMagnelli et all 120091) . We calculate the 
SFR from the UV and IR lumin osities assuming a Kroupa 
IMF, following iBell et alj (120051) . Average SFR values are 
therefore obtained for each subset of galaxies split by red- 
shift in three samples. We combine bootstrap errors with 
scatter between templates to compute errors in estimating 
the SFR. Systematic errors in estimating S FR have n o t been 
incl uded explici t ly, bu t are discussed in IBell et all d2005l) 
and IZheng et al.l ( 120061) . We also compute the average stel- 
lar mass for each subset of sample galaxies using stellar 
mass estima tes obtained from COMBO- 17 optical SEDs by 
iBorch etafl (120061) . Table [2] lists the results, including ob- 
ject number, average SFR, average stellar mass and volume- 
averaged SFR for our samples. 

12 Adoption of the u pdated IR SED templates based on Spitzer observa- 
tions iRieke et al. 2009) gives consistent results within the uncertainties. 
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We further calculate the SFR for sample galaxies that are 
individually detected at 24 fim. The IR luminosities of these 
galaxies typically dominates (> 80%) the bolometric lumi- 
nosity (UV+IR). The 5 cr detection limit of our 24 //m imag- 
ing is 83 //Jy. This limit corresponds to an IR luminosity ~ 
3(10) x 10 10 L o and a SFR roughly ~ 3(10) M yr' atz = 0.6 
(1). Our sample of massive galaxies with SFR > 3 M yr~' is 
therefore incomplete at z > 0.6; a cut of SFR > 10M o yr _I 
is complete up to z = 1 . We bear this in mind and discuss its 
effects on the relevant conclusions. 

We could simply use the SFR density determined from the 
ECDFS alone in what follows, by simply dividing the total 
amount of SFR in these subsamples by the survey volume. 
That approach would give very similar results to those that 
we present later (with the exception of the 0.6 < z < 0.8 bin, 
where a somewhat higher SFR density would be determined 
owing to the overdensity in the ECDFS at 0.6 < z < 0.8), 
and our conclusions would be unchanged. Yet, we choose 
to use a slightly more complex approach that attempts to re- 
duce field-to-field variance by comparing the ECDFS to the 
rest of the COMBO- 17 survey. We adjust the SFR densi- 
ties for the massive galaxy samples by multiplying by the 
ratio of the average stellar mass density derived for all of 
COMBO- 17 to the stellar mass density in massive galaxies 
in the ECDFS: i.e., S FR^ ~ SFR ECD fs X P*,aii/p*,ECDFS- 
Similarly, we adjust the inferred SFR density in spheroid- 
dominated galaxies usin g the stella r mass density in red- 
sequence galaxies fr om iBorch etafl ((2006) as a guide (see 
iMcIntosh et al.ll2005l for the comparison between mass func- 
tions of morphology-selected and color-selected early-type 

galaxies): S f7?„>2.5,all ~ SF/? n> 2.53CDFS x P^red,all/p*,red,ECDFS- 

4.1.4. Local comparison sample 

We augment the intermediate-redshift sample with a sample 
of 2177 local galaxies collected from the NASA/IPAC Extra- 
galactic Database (NED) to assess the corresponding volume- 
averaged SFR at z ~ for the three sub-populations. The 
sample galaxies are selected with a 2MASS A'-band magni- 
tude cut K < 12 in the volume of 1500 < cz (kms^ 1 ) < 3000 
and Galactic latitude b > 30°. About 61% of the selected 
galaxies have redshifts from NED. We believe that the red- 
shift identification is not significantly biased and the sample 
is representative of local galaxies . Details about the sam- 
ple completeness can be found in lBell et all (12005b . Stellar 
masses were estimated from the /^-band absolute magnitude 
assuming a K-bwA stellar mass-to -light ratio of 0.6 M /L 
and a Kroupa IMF (Bel l et al.1 12003b . Of the 2177 sample 
galaxies, 1089 have IRAS 60 and 100 //m detections. The 
IRAS detection limit of ~ 3 x 10~ n ergs cirT 2 s is applied to 
the remaining 1088 galaxies. The total IR luminosity derived 
fro m IRAS ob servations is used to estimate the SFR follow- 
ing iBelll (120031) . The typical error is ~ 0.3 dex for both stellar 
mass and SFR. We calculated the volume-averaged SFR for 
the mass limited subsample, and for the mass and SFR limited 
subsample. The morphological Sersic index parameter is not 
available for this local sample, so the volume-averaged SFR 
at z — is missing for the mass and morphology limited sub- 
sample. We note that the local comparison sample provides 
roughly consiste nt results with the optical ly-selected sample 
from the SDSS dSchiminovich et al.ll2007l) . although a small 
offset exists between the two mainly due to the difference in 
SFR estimator. 

4.1.5. Results 



TABLE 3 

Best-fit parameters for the increases of volume-averaged 
SFR/BHAR of galaxy/AGN sub-populations with lookback time. 



population 


A" 


B a 


Ref. b 


all galaxies 


0.11 ±0.01 


-2.03 ± 0.03 


this work 


massive c 


0.12 ±0.05 


-2.43 ± 0.26 


this work 


massive & n > 2.5 


0.18 ±0.11 


-3.40 ± 0.61 


this work 


massive & 5 FR >3 


0.14 ±0.06 


-2.63 ± 0.28 


this work 


massive &SFR>i0 


0.25 ± 0.06 


-3.46 ±0.31 


this work 


AGN (log Lb i >43) 


0.17 


-2.24 


H07,M04 


AGN (log Lboi >44.9) 


0.21 


-2.73 


H07,M04 


AGN (log Lboi >45.5) 


0.23 


-3.03 


H07,M04 



a Best fit parameters for logpsFR = A?lb +B (galaxies) or logpBH + 
3.11 = Af L B + B (AGN) over the redshift range < z < 1, where 
f L B is the lookback time in units of Gyr. The constant 3.11 is the 
scaling factor 1300 in logarithm. 

b Refe rences — B07: IBell et alj (20071) : LF05: ILe Floc'h et alj 
12005); H07: IHopkins et alj <2007bl) ; M04: IMarconi et all ilOO^K 
c Refer to M, > 2 X 10 10 M s . 

Figure |4] shows the volume-averaged SFR as a function of 
the lookback time for three galaxy subsamples. For compari- 
son, we al so sh ow the volume-averaged SFR density for all 
galaxies ( 92.11 thick gray line). All three galaxy subsam- 
ples have a higher average SFR at larger redshift, consistent 
with the overall galaxy population. A straight line is fit to the 
data points, accounting for the error bars on both axes. 13 The 
line is described by logpsFR = Atis — B, where psfr is the 
volume-averaged SFR and fLB is the lookback time in units 
of Gyr. The best-fit parameters given by the error-weighted 
least-squares fit are listed in table [3] 

Integration of the volume-averaged SFR over < z < 1 
shows that massive galaxies contribute 45% to the overall 
volume-averaged SFR. Most of this SFR (>85%) is con- 
tained in the mass-limited subsample defined to have SFR 
> 3M Q yr _1 ; this reflects that the t ypical SFR of a m assive 
star-forming galaxy is > 3M yr~' (Noeskeetal. 2007|). A 
subsample limited at > lOM yr~' contains 43% of the SFR 
in massive galaxies. 

Star formation in the host galaxies of X-ray-detected AGN 
is not counted. The difficulty in estimating star formation in 
these galaxies is how to separate the IR emission powered by 
star formation from that by AGN. We note that the removal 
of X-ray detected sources in our sample selection does not 
influence our results. Indeed, the 24ytzm fluxes from the X- 
ray detected sources are negligible compared with the total 
24 jum fluxes from thos e X-ray undetected ones (< 15%; see 
also iBrand et al. 2006). The average SFR is dominated by 
galaxies without significant AGN, indicating that SF and BH 
accretion take place in different phases or over different time 
intervals. 

4.1.6. Hosts of star formation: spheroids or disks? 

About one quarter of the SF in massive galaxies is con- 
tained in 'spheroid-dominated' galaxies with n > 2.5; given 
that 45% of the total SFR density at z < 1 is contained in 
the massive galaxy sample, this result is consistent with ear- 
lier results that fo und that < 10% of all SF is in lenticular or 
elliptical galaxies dBell et alj|2005l) . 

Our SFR estimates measure the SF in the whole galaxy. 
The IR data we use do not have sufficient spatial resolution 

13 The vertical error shows the random uncertainties in SFR density. The 
horizontal "errorbar" indicates the redshift range and the data points are uni- 
formly weighted in redshift space. 
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Fig. 4. — The volume-averaged SFR as a function of lookback time for massive (M„ > 2 X 10 10 Mq) galaxies (top-left), massive high-SFR (SFR > SMoyr" 1 ) 
galaxies (top-right) and massive spheroid-dominated (Sersic index n > 2.5) galaxies (bottom-left). Horizontal errorbars represent the redshift range and vertical 
eiTorbars represent 1 o~ errors derived from bootstrapping. The dashed lines show least-squares fits to data points (taking into account the errorbars on both axes). 
The gray thick line in each panel is the least-squares fit to all data points at z < 1 listed in Figure[T] giving best estimates of the global volume-averaged SFR. The 
dotted lines are the same in all three panels, showing the volume-averaged BHAR attributed to luminous AGN (Lb i > 10 44 9 ergs~') converted from the best-fit 
LDDE model of the bolometric LFs given in HRH07 with a radiation efficiency e = 0.1, scaled by a factor of /b e -bh x /recycle ■ Here /bg-bh is the local Mbh /M t 
ratio of 650 given by the local BH-bulge mass relation (Haling & Rix 2004), and / reC ycie is the recycling factor. We adopt / re cycle = 2. Assuming an average 
Eddington ratio log (Lbol/^Edd) = -0.60±0.3 for luminous AGN i Kollmeier et al. 2006), the bolometric luminosity Lbol > 10 449 ergs~' corresponds to BH mass 
Mbh S: 2.6 X 10 7 M Q and accordingly spheroid mass 2 X 10 10 M Q . 



to distinguish between SF in spheroidal and disk components 
within a galaxy. Accordingly, we adopted a simplistic ap- 
proach of assigning all SF in concentrated n > 2.5 spheroid- 
dominated galaxies to spheroids (an overestimate), and as- 
signing all SF in n < 2.5 disk-dominated systems to disks 
(some of that SF may be in spheroids in these disk-dominated 
galaxies; this would drive one towards underestimating the 
SFR in disks). While mid-IR imaging data from JWST will 
help to address this problem in the long-term, our present 
approach is essentially all that the current data support, and 
provides a reasonable guide to the possible contribution of 
spheroids to the volume-averaged SFR density. Bearing in 
mind the caveats, we conclude nonetheless that the bulk of 
the overall SF is associated with disks, rather than spheroids; 
we estimate that < 30% of SF in massive galaxies can be asso- 
ciated with n > 2.5 spheroid-dominated galaxies at any epoch 
in the last 8 Gyr. 



4.2. BH accretion in Luminous AGN 

It would be ideal to measure the BH accretion rate for our 
sample of galaxies in the ECDFS. Unfortunately, episodes of 
intense accretion by supermassive black holes are rare, and 
the number of AGN in the ECDFS is too low to allow an anal- 
ysis as detailed as our above discussion of the SFR density. 
Accordingly, we here adopt an intermediate approach based 
on luminosity-selected samples (derived from larger-area sur- 
veys) that gives a good estimate of the BHAR density in the 
massive galaxy population. 

It is generally assumed that BH accretion obeys the Edding- 
ton limit. A luminous AGN is normally caused by intense 
accretion onto a massive BH, while a faint AGN can be due 
to either a massive BH with low accreti on rate or a le s s mas- 
sive BH with higher accretion rate. Kollm eier et alj (2006) 
found that luminous AGN have an average Eddington ratio 
iboi/^edd - 0.25 with a scatter of 0.3 dex over a wide redshift 
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range. This implies a strong correlation between AGN lumi- 
nosity and BH mass. 

We want to estimate the BHAR in massive BHs in galaxy 
populations of mass M» > 2 x 10 10 M o , and compare it to 
the SFR in the same galaxi es. According to the local BH- 
bulge mass relation given by Hari ng & Rixl (2004), spheroids 
with mass M« > 2 x 10 10 M o host SMBHs with mass M BH > 
2.6 x 10 7 M Q . C ombined with the local mass function from 
iBell et al] (2003), we assume a mean spheroid-BH mass ra- 
tio M,/Mbh = 650 for local massive spheroid-dominated 
galaxies. Assuming Eddington ratio Lboi/^edd = 0.25 and 
radiation efficiency e - 0.1, we infer that massive BHs of 
M B h > 2.6 x 10 7 M power luminous AGN of log L bo i > 44.9 
with BHAR > O.13M yr~'. We have seen that such intense 
BH accretion events statistically match starburst events with 
SFR > 3M yr~', according to the universal agreement be- 
tween BH accretion events and SF events (a factor of 20 in 
intensity; see $3). 

If the local BH-bulge mass relation holds at all cosmic 
epochs, then the mass growth of SMBHs and that of spheroids 

should follow the same relation as = 650. We use 

<M BH > 

the best-fit LDDE model of bolometric LFs from HRH07 to 
calculate the volume-averaged BHAR contained in luminous 
AGN of logLboi > 44.9 over < z < 1. The results are 
scaled up by a factor of 1300 and shown with the dotted 
lines in Figure [4] (same in all four panels). The scaling fac- 
tor 1300 accounts for the mean local spheroid-BH mass ratio 
650 and the recycling factor of two. Note that some studies 
(e.g. iTreu et al] l2007) have claimed to find evidence that the 
BH-spheroid mass ratio was larger in the past (although see 
lLauer et aTl (l2007bl) ). If this is the case, it would obviously 
invalidate the "strong" co-evolution picture. It would impact 
our calculation by changing the slope of the BHAR function 
when cut by the "matching" BH mass, because a given bulge 
mass would correspond to a larger BH mass and therefore to 
a higher AGN luminosity at high redshift. We discuss the im- 
plications of such an effect in Section|5] 

The open circles mark redshifts from z = to z = 1 with 
a step size Az = 0.1. We fit a straight line to these points 
and present the best-fit parameters given by the least-squares 
fit in Table [3] In addition, we also fit the overall BHAR 
(in AGN of log Lboi > 43) and tabulate the corresponding 
best-fit parameters in Table [3] We estimate that luminous 
(log L b oi > 44.9) and very luminous (logL bo i > 45.5) AGN 
are responsible for ~ 70% and ~ 45% of the overall BH ac- 
cretion over 0< z <1, respectively. We caution that the cor- 
respondence between AGN luminosity and BH mass depends 
on the Eddingt on ratio. The uncerta inty of 0.3 dex in the Ed- 
dington ratio (Kollm eier et al.ll2006l) causes an negligible er- 
ror in the sum of the intense BH accretion. The scatter in 
the local BH-bulge relation (~ 0.3 dex) introduces a compa- 
rable error. It is worthwhile to note that the Eddington ra- 
tio may be luminosity dependent: it is possibly close to unity 
for luminous quas ars and increasingly sub-Eddington f or low- 
luminosity AGN dBabic et al. 2007t iBundy et al.ll2008|) . If so, 
our calculations will underestimate the accretion for massive 
BHs. However, as long as the average accretion rate for the 
objects we include in our calculation is similar to that of the 
Kollmeier et al. sample, our results should be reasonably ac- 
curate. In any case, the overall BH accretion is almost cer- 
tainly dominated by luminous AGN associated with massive 
BHs. We emphasize that our conclusions essentially rely on 
the relative density, i.e., the BHAR density of luminous AGN 



to the total, and are marginally affected by how the numbers 
are calibrated. 

4.3. Host galaxy connections 

We showed in 94.1.5l that ~45% of SF happe ns in massive 
galaxies (M, > 2 x 10 10 M o ). We estimated in 94~2l the total 
amount of BH accretion in massive galaxies using an accre- 
tion rate-limited sample, in conjunction with an assumption 
about the distribution of Eddington ratios, to conclude that 
~ 70% of BH accretion should be contained in the massive 
gal axy p opulation. The uncertainties inherent in our analysis 
in 94.21 are considerable, but it appears that the bulk of BH 
accretion tends to occur in massive systems, whereas much of 
the SF happens in lower mass s ystems . 

Furthermore, we showed in 94. 1.61 that < 30% of the star 
formation in massive galaxies (M„ > 2 x 1O IO M ) happens 
in spheroid-dominated systems: star formation, to first or- 
der, happens in disks. Although a substantial fraction of 
AG N by number reside in late- and inter mediate-type galax- 
ies dPierce et al.ll2007t iGabor et al]l2009l) . the luminous ones 
tend to reside in massive, early-type galaxies. Therefore, it 
is clear that the bulk of the black hole accretion happens in 
spheroid-dominated galaxies: it appears that BH accretion is 
happening in different objects or at a different phase in the 
life-cycle of galaxies than the bulk of the star formation. 

We illustrate this point with an example from the 
ECDFS (this exercise simply reinforces the conclusions 
of, e.g.. iGrogin et al] 120051: iNandra et al] 120071: iPierce et al] 
1 20071: iGeorgakakis et al] 120081: lAlonso-Herrero et al] 120081: 
IGabor et al j 120091) . We select massive galaxies with M„ > 
2.5 X IO^Mq in the redshift range 0.4 < z < 0.8 in the 
ECDFS to create a parent sample. From the sample, we se- 
lect those detected i n the 250 ks Chandra X-ray observation 
dLehmer et al] 12003b to make a subsample of galaxies hav- 
ing luminous AGN with soft X-ray luminosity > 10 42 erg s _1 . 
The IR-detected galaxies, i.e. those detected at 24yum, are 
selec ted to comprise a subsample of star-forming galaxies 
(see 94_T3|f or the d etails of the MIPS 24 fim observation and 
Donle y et al] 120081 for further discussion). This subsample 
is also limited to have soft X-ray luminosities < 10 42 ergs~' 
to avoid AGN contamination to the IR. Sersic indices are de- 
rived for all sample galaxies from their GEMS/ACS images 
(HauBleretal. 2007). The X -ray-detected AGN are generally 
obscured in the optical and contribute marginal contamination 
(<~ 10% of central flux) to their host galaxies. Therefore the 
Sersic index is derived from the GEMS/ACS imaging of the 
host galaxies is insensitive to the central AGN. Figure[5]shows 
the histograms of Sersic index n for all galaxies, subsamples 
of IR-detected galaxies and X-ray-detected galaxies. One can 
see that two thirds of X-ray-detected galaxies have n > 3, 
while two thirds of IR-detected galaxies have n <3, sug- 
gesting that the X-ray-detected AGN host galaxies are more 
concentrated than the IR-detected star-forming galaxies. A 
Kolmogorov-Smirnov (K-S) test indicates that the probability 
for both distributions to be drawn from the same distribution 
is less than 10~ 9 , implying that the IR-detected galaxies and 
X-ray-detected galaxies are different populations. Put simply, 
star formation and SMBH accretion take place in dramatically 
different systems: the majority of the global star formation 
happens in late-type galaxies (i.e. disks), while the dominant 
phases of BH accretion occur in intermediate- and early-type 
galaxies. 

5. SUMMARY AND DISCUSSION 



Constraints on the co-evolution of BH and galaxies 



11 



0.3 



0.2 



c 
o 



o 
o 



0.1 



0.0 



— , — i — i — | — i — , — 

all galaxies 

IR detected 

— " X-ray detected 



I 



4 
n 



Fig. 5. — Histograms of Sersic index n for IR-detected galaxies (the thick- 
solid line) and X-ray-detected galaxies (the thick-dashed line) for galaxies 
with M, > 2x 10 10 M and 0.4 < z < 0.8. Here the IR-detected galaxies refer 
to star-forming galaxies detected at 24 /jm and with soft X-ray luminosity < 
10 42 ergs-'. TheX -ray-detected galaxies are those detected in the Chandra 
soft band with X-ray luminosity > 10 42 erg s , suggestive of AGN activity. 
The X-ray-detected AGN host galaxies are more concentrated than the IR- 
detected star-forming galaxies. 

We have shown that the global SF history is proportional to 
the global BH accretion up to (at least) z ~ 1, and the propor- 
tionality factor (2000) is that expected if the BHAR is a fixed 
fraction of the star formation that ends up in galactic bulges. 
This parallel between the cosmic SF history and the cosmic 
BH accretion history is an important constraint on the co- 
evolution between SMBHs and galaxies. It suggests that over- 
all accretion-driven BH growth and overall SF-driven galaxy 
growth trace each other, and that the globally averaged BH 
mass to stellar mass ratio remains constant. Furthermore, we 
find that the SFR and BHAR multiplicity functions are also 
scaled versions of one another at z < 1, suggesting that the 
intensity of star formation and BH accretion are at least sta- 
tistically linked. 

We then refined our analysis to test the hypothesis that star 
formation and BH accretion do not just trace each other in a 
statistical fashion, but that the SFR and BHAR are related by 
a fixed factor f co = /recycle /bg-BH - 1300 in every bulge/BH- 
building event. We focused our analysis on massive galax- 
ies with M» > 2 x 1O 1O M , for which our observational 
sample is complete to z ~ 1. Assuming that the local rela- 
tionship between BH mass and bulge mass remains constant, 
spheroids with M„ > 2 x 1O 1() M should host BH with mass 
Mbh > 2.6 x 10 7 M o . Assuming an ave rage Eddington ratio 
iboi/^Edd = 0.25 dKollmeier et al.| [2006). and e = 0.1, these 
BH, when active, power AGN with BHAR > O.13M yr -1 
or bolometric luminosity logLboi ^ 44.9. We computed 
the BHAR contributed by AGN with bolometric luminosities 
above this limit, and compared it with the SFR contributed 
by galaxies with M t > 2 x 1O 1O M . We found that the 
SFR and BHAR again trace each other but are offset by ap- 
proximately the factor of f m =* 1300 expected for massive 
spheroid-dominated galaxies. This result is nearly unchanged 



if we compute the SFR from galaxies limited in both stellar 
mass and SFR (SFR> 3M yr _1 ), as these "high-intensity" 
SF events dominate the SFR budget in massive galaxies. The 
massive galaxies are simply dominating the total SFR, and the 
big accretion systems are dominating the accretion history (al- 
though with a little more evolution). Just as before when we 
compared the overall SFRD with the overall BHAR we got a 
match, we should get a match for the top two panels. 

Further examination of the nature of the galaxies hosting 
star formation and BH activity showed that only about 30 % 
of the SF "budget" needed to match the BHAR is observed in 
massive spheroid dominated galaxies. Furthermore, the distri- 
bution of Sersic indices for star forming galaxies and that for 
AGN detected in the soft X-ray indicate that these two classes 
of activity tend not to reside in the same types of host galax- 
ies. The bulk of star formation occurs in isolated disks while 
the majority of luminous AGN hosts are intermediate- and 
early- t ype galaxies (see a lso Grogi n et al.ll2005t iNandra et al.1 
20071; iGeor gakakis et al. 2008; Alons o-Herrero et al.l 120081 
Watson, et al. 2009; Gabor et al. 2009). These results clearly 



rule out the simplest picture of object-by-object co-evolution, 
in which BH and bulges grow proportionally and simultane- 
ously. This is also supported by the detection of a small frac- 
tion of luminous X-ray detected AGN in spirals, indicating 
that B H accretion may happe n in non-bulge-building events 
(e.g.. lGeorgakakis et al .1120091) . 

There are both observational claims and theoretical expec- 
tations that the relationship between BH mass and spheroid 
mass might have been larger in the past. While this evolution 
would impact our calculation by shifting the mass of the BH 
(and therefore the cutoff luminosity for AGN) to include in 
our "matched" BHAR, this would only cause a minor change 
in the slope of the BHAR and cannot remove the mismatch 
between SF and BHA in spheroid dominated galaxies. There- 
fore, this would not change any of our conclusions. 

There are two possible ways to reconcile this apparent para- 
dox. One possibility (which shall now call Merger Scenario 1) 
is that bulges and BH grow in proportion to one another, and 
in the same events, but with an offset in time. For example, in 
simulations of galaxy-galaxy mergers including BH growth, 
the luminous AGN phase is delayed by about a dynamical 
time (a few hundred Myr) relative to the peak of the starb urst 
activity dDi Matteo et al.l20"oHlHopkins et al.l2005Ll2006b . In 
this picture, we would expect that most AGN hosts would re- 
semble late-stage mergers, i.e. objects with spheroidal mor- 
phology but young stellar populations, in plausible agreement 
with our observational results. Moreover, in this picture, we 
would expect that a small fraction of AGN to be detected in 
early stage mergers, and the star formation associated with 
this BH growth should be observed in early-stage mergers, 
i.e. close pairs or highly morphologically disturbed objects. 

However, other studies have found that only a small frac- 
tion of th e global sta r formation at z < 1 is o ccurring in galaxy 
mergers (iBell et al.ll200H IWolf et al.ll200l Uogee et al.ll2009[ 
Robai naetal.ll2009F the strongest limit is set by Robaina 
et al., who find that <10% of the star formation in massive 
galaxies (m« > 1O 1() M ) is directly triggered by all identi- 
fiable phases (from close pairs to morphologically disturbed 
remnants) of major galaxy merging. This is strong evidence 
against the Merger Scenario 1 . Therefore, while a small frac- 
tion of the star formation associated with the BH growth that 
we witness at z < 1 may be occurring in the same merger 
events that feed the BH, it seems that this does not represent 
the majority of the SF activity needed to "match" the observed 
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BH accretion. 

The second possibility (Merger Scenario 2) is that most star 
formation takes place in isolated disks. When these disks 
merge, the pre-existing stars are scrambled into a dynamically 
hot spheroidal remnant 14 . The presence of even a modest 
amount of gas leads to the formation of a remnant with higher 
phase -space density and h e nce a deeper gravitational potential 
well dDekel & Coxll2006t iRobertson et afll2006bt ICox et alj 
2006). If the depth of the potential well in the central parts of 
the galaxy determines how large the BH can grow, as in the 
picture of self-regulated BH growth outlined by Hopkins et al. 
(2007a), then the post-merger BH will end up on the BH 
mass-bulge mass relation, despite the fact that the SFR and 
BHAR during the merger event may not necessarily obey the 
universal proportionality factor. 

One can easily see that in this scenario, in many cases the 
BH accretion and the star formation will not trace each other 
during the merger event. For example, in a major merger of 
two pure disk progenitors, the scaled BHAR will be larger 
than the SFR associated with bulge growth. Most of the bulge 
mass will be contributed by the pre-existing stars in the pro- 
genitor disks, and the BH will grow to "catch up" with the 
bulge. In the opposite extreme, in a merger of two galaxies 
with pre-existing massive BH, there may be very little BH ac- 
cretion because the AGN will almost immediately "shut itself 
off". In this case, the SF in the merger might be larger than the 
scaled BH accretion (if the progenitors contain enough gas). 

Because what we see is more BH accretion than we would 
expect based on the amount of star formation in spheroids, we 
conclude that most of the mass in these spheroids was con- 
tributed by the pre-existing stars in the progenitors rather than 
by new stars formed in the merger-associated bursts. Indeed, 
we have independent lines of evidence that this must be the 
case. We know that the colors of elliptical galaxies and bulges 
at z < 1 are characteristic of fairl y old stellar populations 
(Galla zzi et al.1 120061 lRenzinll2006L and references therein), 
and similarly, fossil evidence from line strengths in nearby el- 
lipticals also precludes the very recent forma t ion of a large 
fracti on of the stellar mass (Trager et al. 2000; Thom as"et alj 
120051) . In addition, this picture is consistent with the rela- 
tively small fraction of star formation associated with ongo- 
ing mergers at z < 1 : in cosmological mod els that incorpo- 
rate a scenario like Merger Scenario 2 (e.g. ISomerville et alj 
2008), the predicted fraction of star formation triggered by 
major mergers at z < 1 is about 7%, in excellent agreement 
with the observational estimates (see Robaina et al. 2009 for 
details of the comparison). 

A slightly modified version of Merger Scenario 2 can 
also be accommodated in a Universe in which the BH- 
spheroid mass ratio has decreased with time. This is in 
fact expected in the theoretical picture that has emerged 
from the analysis of hydrodynamic simula tions of galaxy 
mergers with self-regulated B H growth (e.g. Robertso n et alJ 
l2006at iHop kins et al. 2007a). These studies find that the 

14 In what follows, much of what we argue for merger scenario 2 would 
also apply if the main route for forming low-mass spheroids were secular 
evolution. The key point for this paper is that the bulk of the stars that end 
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final BH-spheroid mass ratio depends on the gas frac- 
tion of the progenitors. This is because gas-rich progen- 
itors can dissipate energy during the merger, and all else 
equal, form m ore c ompact remnants than gas-poor mergers 
dDekel & Coxll2 006; Rob ertson et ai1l2006bt :lcbx et alj|2006l; 
ICovington et al.l 120081) . In this deep potential well, the BH 
can grow to a larger mass (relative to the spheroid) before the 
pressure-driven outflow halts further accretion. If gas frac- 
tions in merging galaxies were higher at high redshift, as is 
generally expected, this would then lead to the prediction that 
BH were more massive relative to their spheroids at high red- 
shift. A significant fraction of mergers at z < 1 are expected 
to involve at least one galaxy that already contains a massive 
BH. If the pre-existing BH is already larger than the "criti- 
cal mass" set by the new remnant, further BH growth will be 
stifled. Theoretical models predict that this expected evolu- 
tion is relatively mild : less than a factor of two since z ~ 1 
dHopkins et al.l | 2007al) . 

Bor ch et all ([2006) have shown that the stellar mass in red 
(predominantly spheroidal) galaxies has increased by a fac- 
tor of 2-3 since z ~ 1, while the stellar mass in blue (disk- 
dominated) galaxies has remained about the same over this 
time period. At the same time, we know that the bulk of new 
star formation is occurring in blue (disk-dominated) galax- 
ies, implying that mass in blue (disk) galaxies must be trans- 
formed into red (spheroidal) galaxi es at approxima tely the 
same rate that new stars are forming (Bell et al. 2007). Given 
that most of the red galaxies are s pheroid-dominate d, which is 
well established at least at z < 1 dBell et al.ll2004al) . the mass 
on the red sequence cannot grow simply via the quenching 
of star formation in disks — a strong dynamical process like 
merging is required. 

Taken together, this evidence strongly favors a picture in 
which, at least at z < 1, the bulk of the growth of mass in 
spheroids is due to merging of pre-existing, already old disk 
stars, and BHs preferentially grow to "catch up" to their newly 
assembled bulges following a merger event. In this case, the 
parallel evolution of the global SFR and BHAR may simply 
reflect the apparent coincidence that, at z < 1, the rate of 
formation of new stars is approximately equal to the rate that 
stellar mass is transferred from disks into bulges via mergers. 
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up in the spheroid are formed long before the epoch of actual bulge creation, 
which applies equally in Merger Scenario 2 and secular evolution. 
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